AnĂĄlise para varias variĂĄveis

Entrada de dados

library(tidyverse)
data_set <- read_rds("../data/aula6-01.rds")
glimpse(data_set)
## Rows: 45
## Columns: 10
## $ macrofita <chr> "BRASU", "BRASU", "BRASU", "BRASU", "BRASU", "BRASU", "BRASU

## $ dose_t_ha <dbl> 13.5, 13.5, 13.5, 13.5, 13.5, 27.0, 27.0, 27.0, 27.0, 27.0, 

## $ repeticao <dbl> 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 

## $ mo        <dbl> 8, 7, 8, 7, 7, 10, 10, 9, 9, 9, 9, 10, 10, 8, 10, 8, 8, 8, 7

## $ p_h       <dbl> 5.8, 5.8, 5.9, 6.0, 6.0, 6.1, 6.0, 5.7, 6.1, 6.1, 6.1, 6.1, 

## $ p         <dbl> 3, 4, 4, 3, 4, 4, 4, 3, 5, 4, 4, 3, 3, 3, 3, 8, 7, 7, 7, 10,

## $ k         <dbl> 1.01, 0.86, 5.64, 4.57, 4.71, 11.86, 11.95, 7.78, 9.33, 7.70

## $ ca        <dbl> 17.95, 19.95, 18.98, 15.81, 16.37, 21.79, 21.03, 17.19, 20.6

## $ mg        <dbl> 6.29, 6.90, 6.64, 5.57, 5.76, 7.50, 7.55, 6.14, 6.91, 5.93, 

## $ al        <dbl> 0.28, 0.25, 0.19, 0.26, 0.20, 0.21, 0.22, 0.22, 0.20, 0.14, 


Tabelas com DT

library(DT)
datatable(data_set)

DICA:PARA COLOCAR NEGRITO, COLOQUE A PALAVRA ENTRE 4 ASTERISCOS

DICA: PARA COLOCAR EM ITÁLICO, COLOQUE A PALAVRA ENTRE DOIS ASTERISCOS

Selecionar os nomes das variĂĄveis

variaveis <- data_set %>% 
  select(mo:al) %>% 
  names()
variaveis
## [1] "mo"  "p_h" "p"   "k"   "ca"  "mg"  "al"

Fazer anĂĄlise para todas as variĂĄveis

for(i in seq_along(variaveis)){
  my_var <- variaveis[i]
  macrofitas <- data_set %>% pull(macrofita)
  doses <- data_set %>% pull(dose_t_ha)
  y <- data_set %>% pull(my_var)
  print("------------------")
  print(my_var)
  print("------------------")
  media <- mean(y,na.rm=TRUE)
  desv_pad <- sd(y,na.rm = TRUE)
  cat("Média:",round(media,2),"\n")
  cat("Desvio-padrĂŁo:",round(desv_pad,2),"\n")
  
  # GrĂĄfico de histograma
  plot_hist <- tibble(macrofitas,doses,y) %>% 
    ggplot(aes(x=y))+
    geom_histogram(bins = 7,fill="grey",colour="black")+
    theme_bw()+
    labs(x=my_var)
  print(plot_hist)
  
  # AnĂĄlise de variĂąncia
  ExpDes.pt::fat2.dic(macrofitas,doses,y,
                      sigF = 0.15,sigT = 0.05,
                      mcomp = "tukey",
                      quali = c(TRUE,FALSE),
                      fac.names = c("macro","dose"))
  
}

Modelo do matemĂĄtico do experimento fatorial em DIC

\[ y_{ij}=\mu + \tau_{i} + \beta_{j} + (\tau \beta_{ij})_{k} + \epsilon_{ijk} \]

Quadro da anĂĄlise de variĂąncia

FV GL
A \((a-1)\)
B \((b-1)\)
AB \((a-1)(b-1)\)
R \(ab(r-1)\)
T \(abr-1\)

Adição da imagem do prédio central

GrĂĄfico interativo

library(plotly)
graf <- data_set %>% 
  group_by(dose_t_ha,macrofita) %>% 
  summarise(mo=mean(mo)) %>% 
  ggplot(aes(x=dose_t_ha,y=mo,fill = macrofita))+
  geom_col(position = "dodge",colour="black")+
  theme_bw()+
  scale_fill_viridis_d()
ggplotly(graf)

Mind blowing (meme)

## Typing cat (meme)

ReferĂȘncia